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We study theoretically resonant inelastic x-ray scattering (RIXS) at the Cu La-edge 
in a typical parent compound of high-Tc cuprate superconductors La 2 Cu 04 on the 
basis of a detailed electronic band structure. We construct a realistic and precise tight- 
binding model by employing the maximally-localized Wannier functions derived from a 
hrst-principles electronic structure calculation, and then take account of the Coulomb 
repulsion between d electrons at each Cu site. The antiferromagnetic ground state is 
described within the Hartree-Fock approximation, and take account of electron cor¬ 
relations in the intermediate states of RIXS within the random-phase approximation 
(RPA). Calculated RIXS spectra agree well with the experimentally observed features 
including low-energy magnon excitation, d-d excitations, and charge-transfer excita¬ 
tions, over a wide excitation-energy range. In particular, we stress the importance of 
photon polarization dependence: the intensity of magnon excitation and the spectral 
structure of d-d excitations depend signihcantly not only on the polarization direction 
of incident incoming photons but also that of outgoing photons. It is demonstrated that 
the single-magnon excitation intensity is maximized when the polarization directions of 
incoming and outgoing photons are perpendicular to each other. 


1. Introduction 

Resonant inelastic x-ray scattering (RIXS) provides a powerful optical method of 
observing elementary excitations in solids. Particularly, recent tremendous progress in 
energy-momentum resolution and theoretical understanding is promoting RIXS to one 
of major techniques of measuring various electronic and magnetic excitations. Among 
RIXS phenomena, RIXS at the transition-metal edges is attracting much interest, which 
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is much suitable for elucidating electronic excitations in strongly correlated transition- 
metal compounds including the high-T^ cuprate superconductors and various magnetic 
materials. It is notable that RIXS can provide us a variety of electronic excitation spec¬ 
tra, depending on which absorption edge of transition metal is utilized. At the A'-edge 
in the hard x-ray regime, an inner-shell Is electron is promoted to conduction p bands. 
Accompanied with this promotion, strongly correlated d electrons near the Fermi level 
are excited to screen the created inner-shell Is hole. Thus correlated d-electrons are indi¬ 
rectly excited. In this case, electronic excitations should conserve the total spin moment 
of the d-electron system, and therefore are restricted only to charge-orbital excitation 
processes or to generation processes of even-number magnons. In fact, charge-transfer 
excitations,^’ d-d orbital excitations,®^ and two-magnon excitations®^ are observed in 
transition-metal compounds through iF-edge RIXS. On the other hand, at the L-edge 
in the soft x-ray regime, an inner-shell 2p electron is promoted to correlated conduction 
d bands. Thus correlated d-electrons are directly excited. Since the spin and orbital 
angular momenta of the 2p states polarize each other due to strong spin-orbit coupling, 
the d electrons can be excited not only in the charge-orbital channel but also in the 
spin channel. In application to high-T^ cuprate superconductors, in which the spin mo¬ 
ment is induced on the dx' 2 -y 2 orbital, it had been believed that spin-flip scattering is 
not allowed in L-edge RIXS.^’®^ Actually, more recently Ament and collaborators ver- 
ihed that spin-flip excitation is indeed allowed when the spin moments point parallel 
to the basal plane.In fact, magnon excitation is clearly observed in copper oxides 
where the spin moments are aligned along the basal plane.Magnon excitation is 
observed persistently even in doped metallic cuprates,and L-edge RIXS promises 
to be a new technique of elucidating the pairing glue in high-Tc superconductivity. The 
observed energy dispersion of the magnon peak agrees well with neutron scattering 
experiments. i4, 15 ) 

In our present work, we discuss theoretically RIXS at the transition-metal L-edge. 
Intensive theoretical studies have been performed so far to analyze the L-edge RIXS 
in strongly correlated transition-metal compounds.®’In most of precedent studies, 
simplihed effective models such as impurity Anderson models, hnite-size cluster mod¬ 
els, or Heisenberg models have been adopted to describe the electronic structure or the 
magnetic ground state, and the energy window and the number of momentum points 
allowed for calculation are inevitably restricted. Heisenberg antiferromagnetic superex¬ 
change couplings are often treated as tunable parameters to £t to experimental data. 


/ 
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and microscopic quantitative grounds for the values of coupling parameters are lacking. 
Therefore, it will be meaningful to develop another approximate but useful calculations 
without the above drawbacks, for analyzing experimental data in detail. 

In the present study, we present a microscopic theoretical formulation of RIXS 
at the L-edge and apply it to a typical parent compound of high-Tc superconductors 
La 2 Cu 04 . To describe the electronic structure of La 2 Cu 04 precisely, we use maximally 
localized Wannier functions (MLWF) derived from hrst-principles electronic structure 
calculation,^^) and determine the antiferromagnetic ground state within the Hartree- 
Fock (HF) approximation. Electron correlations in the intermediate states are treated 
within the random-phase approximation (RPA). Our approach is based on perturbation 
expansion in Coulomb interaction, and a natural extension of our previous theoretical 
formulation of iF-edge RIXS.^^’^®) It is approximate one but applicable to realistic 
and complex electronic structures such as multi-orbital systems. Similar formulation is 
already developed in Ref. 27. 

The article is constructed in the following way: In § 2, we present our Hamiltonian 
and theoretical formulation to calculate RIXS spectra. In § 3, we present calculated 
RIXS spectra of single-magnon and d-d excitations, and their dependences on polar¬ 
ization direction and scattering angles, and compare some of them with experimental 
data. In particular, we shall see the intensity of magnon excitation and the spectral 
structure of d-d excitations depend signihcantly not only on the polarization direction 
of incident photons but also that of outgoing photons. In § 4, some discussions and 
remarks on our formulation and results are given. In § 5, the article is concluded with 
brief comments. 

2. Formulation of RIXS 

2.1 Hamiltonian 

To discuss the RIXS process microscopically, we use the following form of Hamilto¬ 
nian: 

H = Hn.f. + H2p + H2p-d + Hx, ( 1 ) 

where H^j, describes the correlated electrons near the Fermi level, and H 2 p and 
describe the inner-shell 2p electrons and the dipole transition by x-rays, respectively. 
H 2 p-d is the Coulomb interaction between the 2p and transition-metal d electrons (Cu- 
Sd electrons in the case of La 2 Cu 04 ). We present details of each term in the following. 
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To construct the Hamiltonian part Hnj, for La 2 Cu 04 , firstly we perform first- 
principles band structure calculation assuming the paramagnetic state.To express 
the electronic orbital bases and scattering geometry, we take the coordinate system 
where the principal axes of nearest-neighbor Cu-0 bonds are parallel along the carte¬ 
sian axes (see Fig 1(a)). Then we perform tight-binding htting to the obtained energy 
bands near the Fermi level by using the wannier90 code,^'^^ where we take hve d or¬ 
bitals at each Cu, and three p orbitals at each O site. Thus we include 17 MLWF’s 
in the unit cell, since there are one Cu and four O sites in the unit cell. The orbital 
bases are dehned in terms of the coordinate axes in Fig. 1(a), and the spin states are 
specihed with respect to the 2 ;-axis. Thus we obtain a tight-binding model to £t the 17 
bands near the Fermi level (See Appendix A for details). Adding the on-site Coulomb 
interaction part, we have the Hubbard-type Hamiltonian Ffn./. in the following form: 




t.m. @ri 


5EEE^ 






a'^'^^3 ^ii 4 a: 


i uu' 


( 2 ) 

where and are the electron creation and annihilation operators for orbital £ 
with spin a at site i. Throughout the present study, we always work with the electron 
representation not with the hole representation, ‘t.m.’ in the summation with respect to 
i means summing only over transition-metal sites. is fh® on-site 

Coulomb integral at transition-metal (i.e., Cu) sites. In the summation with respect to 
i, at the top means orbital £ should he on the site r^. One-particle energy at 

orbital £ is given by ei = = 0). We modify the one-particle energy ei for Cu-d 

orbitals, to obtain a realistic level scheme of the local Cu-d orbitals and d-d excitation 
energies, as explained in Appendix A. Hereafter we use the following convention for later 
discussions: if £ denotes a d orbital at a Cu site (e.g., £ = xy), then = da^ = di^, if £ 
denotes a p orbital (e.g., £ = x) at an O site, then = pua, and so on. We expect that 
confusion between Cu-p and 0-p orbitals will not occur, since the 0-p orbitals do not 
appear explicitly in the following discussions. Here we introduce the values of on-site 
Coulomb interaction ^^'Ch Cu site in the form of Slater-Condon integrals 

(see Ref. 29 for the dehnition of Slater-Condon integrals and their relation to hx,i 2 -/ 3 M)'- 
F^{d,d) = 6 eV, F‘^{d,d) = 11.5 eV, F'^{d,d) = 7.4 eV. The values of F‘^{d,d) and 
F^(d, d) are taken from Ref. 30, and F^{d, d) is determined to reproduce the insulating 
gap about 2.2 eV.^^^ Our choice of these Coulomb integrals corresponds approximately 
to 7/ ~ 7.5 eV, U' ~ 5-6 eV, and J ~ 0.6-1.2 eV, where U, U' and J are the intra-orbital 
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and inter-orbital Coulomb repulsions and the Hund’s coupling, respectively. For Hnj., 
we determine the antiferromagnetic ground state with the spin moments m || [110] 
within the HF approximation (see Appendix A about details of HF calculation). 

For the 2p electrons, we assume completely localized 2p orbitals at each transition- 
metal site (Cu site in the present case of La 2 Cu 04 ): 


t.m. r@ri @ri 

H2p = EE ^2pi'f’i)PimaPimcr + EE ^2p{,f‘i)lmm' ' ^crcr'PimaPim’a' j (3) 

i mm' acr' 

where = ^ 2 p is the one-particle energy of the 2p state, p\^^ and Pima are the 

creation and annihilation operators of 2p electrons with spin a and angular momentum 
m {m = —1, 0,1) at transition-metal site i, respectively. The second term of the right- 
hand side of eq. (3) is the spin-orbit coupling. are the matrix elements of orbital 
angular momentum, and are related to the Pauli matrices as s = cr 12. H 2 p is 
straightforwardly diagonalized, and the 2p eigenstates are characterized by the total 
angular momentum J and its 2; component M, where J = 1/2, 3/2 and M = — J, — J + 
1,..., J — 1, J. The energy levels consist of two kinds of levels as a result from splitting 
due to the spin-orbit coupling: = ^ 2 p - and e 2 pj^ 3/2 = ^ 2 p + 6p/2. The 

former (latter) is two-fold (four-fold) degenerated. In the present case of Cu, we take 
e 2 p = —936 eV with respect to the Fermi level, and ^ 2 p = 13.3 eV, which causes about 
20 eV energy splitting between the L 2 and L 3 edges. The diagonalized eigenstates are 
expressed as 


J, M) = ^ ^ |m, a)Uma-^jM- 

m=lfl,—l (T=t,t 


Specihcally, for the J = 1/2 doublet, 
|J = 1/2,M = 1/2) : 

|J = 1/2,M =-1/2) : 

and for the J = 3/2 quartet, 

\J = 3/2, M = 3/2) 

\J = 3/2, M= 1/2) 

\J = 3/2, M = -1/2) 

\J = 3/2, M = -3/2) 



(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
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H 2 p-d is given by 

t.m. @ri @ri 

H2p-d V2p-d{ri-, maija2-, i'a^, im' (74,1 ( 11 ) 

i ii' mm' <Ji'^cr 4 

where V 2 p-d{ri] ■ ■ ■) is the on-site Conlomb interaction between the 2p and d orbitals at 
transition-metal site rj. i and denote the d orbitals: i, i' = xy, yz, xz, x"^ — y'^, Sz"^ —r"^. 
For convenience in later discnssions, we dehne the following matrix elements: 


V2p-d{ri;JMX;C, J'M') = ■jM^ 2 p-d{ri] mai,ia; i'a\ m'(j 2 )um'a 2 -,J'M>- (12) 


mm' <7ia2 


For nnmerical calcnlations on La 2 Cn 04 , we take F^{p, d) = F‘^{p, d) = 1 eV at Cn sites, 
bnt this interaction is almost irrelevant to spectral weights for the case of La 2 Cn 04 as 
we see later. 

Hx describes resonant 2p-d dipole transition indnced by x-rays: 

t.m.u ^ @ri 

Hx = EyEEE q, e)a 9 e 4 +q£aPfc ma h.C. , ( 13 ) 

i k,q £,m a 


where N is the nnmber of nnit cells, the snmmation in i with ‘t.m.n’ means it is re¬ 
stricted to transition-metal sites in the nnit cell, d],^^ [pkma] is the creation [annihilation] 
operator of transition-metal de [2pm] electrons in the momentnm representation, and 
age is the annihilation operator of a photon with momentnm q and polarization e. We 
assnme the matrix elements of Wi^rniffi Q, e) are given in the form: 


Wi^mirp q, e) 



{de\e ■ p\2pm) 

- — \l - e2pX{M^ ' '^|2pm) 


—ie\/2TT\q\F'''''''{di\e ■ r\2pm) 


(14) 


in natnral nnits (c = h = 1), rrie and e are the mass and charge of an electron, re¬ 
spectively. Vi is the transition-metal site where the x-rays are absorbed or emitted, and 
we assnme that we^rnifi^qje) does not vanish only when both of the orbitals 2pm and 
di lie on the site r*. {di\e ■ r\2pm) can be calcnlated nsing atomic wave fnnctions. For 
convenience in later discnssions, we dehne 

^,JQi ^) = i(T,JMij"i) Qi ®) = ^ ^ ^i) Qi ^^'^m,cr\JMi 

m=l,0,-l 

where C, is orbital-spin nnihed index: ^ = (•^, a). 


(15) 
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Fig. 1. (Color online) (a) Crystal structure of La 2 Cu 04 and MLWF for the state, (b), (c) 

and (d) Definition of the scattering geometry, q and q' are the momentum vectors of the incoming and 
outgoing photons, respectively. 0, (j) and ^ [9\ </>’ and are the Bragg, azimuthal and polarization 
angles, respectively, for incoming [outgoing] photons. In (c), the angle between transferred momentum 
Q and the 2 -axis is approximately given by i5 ~ (0 ~ ^0/2- In (d), e and e' are the polarization 
vectors of the incoming and outgoing photons, respectively. The polarization angle is measured with 
respect to the scattering plane, i.e., tp = 0 [tp = 7r/2] means that the polarization direction is parallel 
[perpendicular] to the scattering plane. 


2.2 RIXS intensity 

RIXS intensity can be obtained by calculating the number of photons generated in 
different states from the incident-photon state per unit time. To do this, we employ the 
Keldysh perturbation theory as in Ref. 32. The present formulation of L-edge RIXS 
is a straightforward extension of our previous formulation on iC-edge RIXS.^^’^®^ The 
RIXS intensity is generally expressed by the diagram (I) in Fig. 2, if assuming that only 
a single electron-hole pair remains in the hnal state. The analytic expression of RIXS 
intensity is obtained from the diagram (I) of Fig. 2 as: 




fcl aia2 


t.m.u. @ri J J' 

EEXE X 

i CC' JJ' M=-J M'=-J' 
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Fig. 2. (I) RIXS intensity represented within the Keldysh perturbative formulation. The wavy 

lines and shaded rectangular represent the photon propagators and electron scattering vertex func¬ 
tion F{ri\q,q'), respectively. A pair of oriented solid lines represent the off-diagonal elements of 
the Keldysh Green’s function, and connect the upper normally-time-ordered and lower reversely- 
time-ordered branches. (II) Approximate expansion for the scattering vertex function F(ri;ki]q,q'): 
(a) F^'^\ri\ki\q,q') for ‘Oth-order process’ (fluorescence-like), (b) F'^p'> {ri\q,q') for ‘p-process’, (c) 
F^‘^\rf,q,q') for ‘ d-process’. The filled triangle and square are the three-point and four-point vertex 
functions to be renormalized by electron correlations, respectively. In (b), the dashed line represents the 
core-hole potential V 2 p-d- Thick solid lines represent the propagator of the inner-shell 2p electrons. (Ill) 
RPA diagrams for the three-point and four-point vertex functions (Q) and (Q) ((a) 

and (b), respectively), where empty squares represent the antisymmetrized bare Coulomb interaction 
^C 2 ^C 4 C 3 Cl ^ electrons at transition-metal sites. 


. ( 16 ) 


J^M']ai,a2 Q-) Q ^ 
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where (/c) is the Keldysh Green’s function, 01^2 are indices for the diagonal¬ 
ized bands, C = (A o'), C' = Ecc' = and ki = (a;i,fci). q 

and q' are the four-momenta of the incoming and outgoing photons, respectively: 
q = (a;, g), q' = {uj',q'). e and e’ are the unit vectors pointing along the po¬ 
larization direction of the absorbed and emitted photons, respectively. Q is the 
energy and momentum loss of the photon: Q = q — q' = {u — u',q — q') = 
(G,Q). Fc,JM,c,'j'M'-,ai,a 2 iTi')ki]q,q') is the scattering vertex function expressed us¬ 
ing only the ordinary causal electron Green’s functions and electron-electron inter¬ 
action. At this stage, we omit uji dependence of K^jm, c'a 2 (^ij O') ^0; i-®-) 
F(jM,CJ'M'm,a 2 {ri; h; Q, q') = F(;jM,CJ'M'-,ai,a 2 iri; ki; q, q'), because it is justihed within 
the following approximation for K^jM,c'J'M';ai,a 2 (’"b q')- Within the HF approxima¬ 

tion, the Green’s functions G^(/c) are given by 

G+(fci) = 27rina{ki)6{ui-Ea{ki)), (17) 

= -2m[l-naiki)]6{uji-Ea{ki)), (18) 


where Ea{ki) is the energy of diagonalized band a, and Uaiki) is the electron occupation 
density i.e., the Fermi distribution function, at momentum ki in band a: na{ki) = 
l/(exp[i7a(fci)/T] -|- 1). Substituting eqs. (17) and (18) into eq. (16), we have 

Oyr 

^{q,e;q',e') = — (A:i)[l — na2{ki -|- Q)](5(r2 -|- Eaj^{ki) — 77^2(fci + Q)) 

ki a\a2 


X 


t.m.u. ®ri 


J 


EEEE E F)wQi^jij^i{ri] q ,e) 

i CC' JJ' M=-J M>=-J> 

2 


(19) 


For calculation of ^ij O'; use perturbation expansion with re¬ 

spect to electron-electron interactions. There are three possible major contributions to 
F(jM,CJ'M'-,ai,a 2 {'^uki;q,q'). The hrst is the zeroth-order term represented by the dia¬ 
gram (II)-(a) in Fig. 2. This diagram presents the main contribution from fluorescence 
processes, where a valence d electron goes down to the 2p state without interacting 
with other electrons. We refer to this contribution as ‘Oth-order process’. The second 
originates from the screening of the 2p core hole. Within the Born approximation with 
respect to the core-hole potential V 2 p-d, this process is expressed by the diagram (II)- 
(b) in Fig. 2. We refer to this contribution as ‘p-scattering process’ or ‘p-process’. The 
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third describes the scattering processes involving both the excited conduction d elec¬ 
tron and valence d electrons. This contribution is expressed by the diagram (II)-(c) in 
Fig. 2. We refer to this contribution as ‘d-scattering process’ or ‘d-process’. Of course, 
in higher-order contributions, more complex diagrams can appear, which cannot simply 
be classified to ‘p-process’ or ‘d-process’. Nevertheless, this classification turns out to 
be convenient for microscopic analysis of RIXS spectra. Thus, we obtain the following 
approximate expression for the scattering vertex function: 

Fc_JM,Q'J'M'\ai,a2iTi) d) d ) = JM,C' ,02 d) d ) 

C1C2 

+ d) d )]) (20) 


where is the diagonalization matrix of the HF Hamiltonian given by eq. (A-3). 

Cn is orbital-spin unihed index: = [(.ni o'n), and where £n represents 

d orbitals at transition-metal sites. For the antiferromagnetic ordered state, the energy 
bands are doubly folded. Therefore, in addition to Cn, a two-valued index is necessary 
to specify which of doubled bands each state is on. But we suppress it for simplicity of 
notation. Contributions from the above three processes are given by 


P 


( 0 ) 




(r*; fci; g, q') = 6 jj/Smm' 


^ + ^2pj{'f’i) — Ea2 {ki + Q) ’ 


( 21 ) 


@ri 


-^C^A^,C'J'M';Ci,C2(^b QiQ) — y2p-d{'f'i] JM, C,^] (^ 4 , J M ) A^2,C4;C3,Cl (Q) 

C 3 C 4 


X 


X 


a k 

'^df(i),a(^)'*^c«c' h)>a(^) 

[cn + e 2 pj{ri) - Ea{k)][u' + e2pj,{ri) - Ea{k)] 


,( 22 ) 


-^CJM,C'J'M';Ci,C2(^d d, d ) — F^2,C4;C3,Cl (Q) 'l^asik + Q)] 


C 3 C 4 


N 


0304 


X- 


'^dQ(i),aS^ Q)'^C3,a3(^ + Q)^C4,a4 (^)^'^C'h)>“4 (^) 


X 


^ + ^2pj (n i) — Ea^ {k + Q) 

l-na^{k) na^{k) 


+ ^2pj,{'^i) ~ Ea^{k) n -|- Ea^{k) — Ea^{k Q) + ie J 


(23) 
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where A^ 2 ,C 4 ;C 3 >Ci (Q) are the three-point and four-point vertex func¬ 

tions, which are represented by the hlled triangle and square in Fig. 2 (II) (b) and (c), 
respectively. Index = di„{i) denotes the d^ state at transition-metal site r* with 
spin a. e 2 pj{ri) = + *r 2 p, where r 2 p is the damping rate of the 2p core hole 

and set to 0.3 eV in the present study. Summations in i with ‘t.m.u.’ at the top means 
that Vi should be restricted only to transition-metal sites in the unit cell, e in eq. (23) 
ensures the causality, i.e., that an electron-hole pair is created only after the resonant 
transition of a 2p electron to d states, e has the physical meaning of the damping of the 
excited electron-hole pair. In our numerical calculations, we set e = 20 meV. 

The vertex functions introduced above are renormalized by electron correlations. 
We take account of electron correlations within RPA. RPA for A^^ (Q) and 

P^ 2 ,C 4 ;C 3 ,Ci(Q) represented diagrammatically in Fig. 2 (III) (a) and (b), respectively. 
The analytic expressions for these diagrams are 


Ac.,C4;C3.Ci(Q) = '^C2C3'^C4Cl 


A//*/ A/ A/ 


rc,.c4;C3,c.(Q) = (25) 

cia (3<i 

where is the antisymmetrized bare Coulomb interaction given by = 

-^q/2;4y4'^<Tia4<^f72<73 - hui 2 M 3 ^<^i<T 3 ^<T 2 <Ti- x{Q) IS the polarization function calculated by 


XC3,C2;Ci,t4(Q) “Cl,a(^)“C4,a(fc)^<3,a'(fc + Q)ul^,a'ik Q)Xa,a'{k- Q), 


Xa,a' {k, Q) 


Ua'ik + Q) - na{k) 


(26) 

(27) 


n -I- Ea{k) — Ea'{k + Q) +ie^ 
where e is interpreted as the damping rate of the excited electron-hole pair near the 
Fermi level, as already introduced above. Solving eqs. (24) and (25), we can determine 

Ac2,C4;C3,Ci(< 5) and F^2,C4;C3,Ci(<5) within RPA. 

To resolve into the contribution from each process of the Oth-order, p-scattering and 
d-scattering, we introduce the process-resolved spectra as follows: 

hF(°)(g,e;g',e') = vEE" ai(fcl)[l ~ na^ikl + Q)]d(r2 -|- Ea-^{ki) — Ea2{ki -f- Q)) 


h\ aia2 
t.m.u. @ri 


J J' 

EEEE E 

i C 1 C 2 JJ' M=-J M'=-J' 


X 
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(28) 


2tt 


W^P\q,e-,q',e') = fEE"«. (fci)[l — na2{ki + Q)]5(fi + Eaiiki) — Ea2{ki + Q)) 

hi aia2 

t.m.u. J J' 

EEEE E Qi ,J' ^') 

i CC' JJ' M=-J M'=-J> 

2 


X 


C1C2 


(29) 


27r 


)fEE"« (fci)[l — na2{ki + Q)](5(r2 + Ea^iki) — £’^ 2(^1 + Q)) 


X 


X 


k\ aia2 
t.m.u. @ri 


J J' 


EEEE E Q, Q 1 ^) 

i CC' JJ' M=-J M'=-J' 

2 

,C'J'M';Ci,C2(^*’ ) 

C1C2 


(30) 


These are obtained from eq. (19) by keeping only one of fep g, g'), 

-^CJM,c'j'M';Ci,C 2 (’’b Q, q') and -F^jLf,c'J'M';Ci,C 2 (’’b Q, q') and setting the rest two to zero in 
eq. (20). 

Further to resolve orbital-excitation processes involved in the Oth-order, p-scattering 
and d-scattering processes, we introduce the orbital-spin-resolved spectra as follows: 


Wl^\^^{q,e]q',e') = 


2-k 


(fci)[l — na2{ki -|- Q)](5(f2 -|- Ea-i^{ki) — Ea2(ki + Q)) 

ki 0102 

J 

EE E 


X 


JJ' M=-J M>=-J> 

m 


^-^C2-/M,CiJ'M';ai,02(^h ^Ij d) d ) 


(31) 


»(‘rE(9,e;9',e') = ^ E E ('“i + <3)^^ + (*=i) “+ <3)) 


X 


e; q\ e') = 9- E E - n«(fci + < 5 )li 5 (n + E„(k,) - E„(k^ + Q)) 


N 

hi aia2 

t.m.u. J J' 

EEEE E Q, Q'y ^') 

i CC' JJ' M=-J M'=~J' 

2 

^^C 2 ,02(^1 Q)'^Ci.«i (^i)-^cjM ,C'J'M';Ci,C 2 (’’b d) d ) 

27r 


(32) 


N 


hi aia2 
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t.m.u. @ri J J' 

xEEEE E q, e)w*(-, J!]^,{ri] q', e') 

i CC' JJ' M=-J M'=-J' 

2 

^“C2,a2(^l + Q)“Cl,ai(^l)^CJM 


(33) 


where, in every of eqs. (31), (32) and (33), the orbital-spin index (^i [^ 2 ] represents the 
orbital and spin of the electron removed below Ep [left above Ep] in the hnal state at 
t = -|-cx), where Ep is the Fermi energy and t is time. Throughout the study, we specify 
spin states by taking the spin axis parallel to the 2 ;-axis. Equations, (31), (32) and (33) 
are obtained by suspending the summation with respect to orbital-spin indices and 
(2 in the right-hand side of eqs. (28), (29) and (30). In eq. (31), if we select orbitals 
Cl and C 2 specihcally, then atomic site r* in the unit cell is determined uniquely to give 
non-vanishing contributions, since the factor does not 

vanish only in the case where both of Ci and C 2 reside on atomic site r*. Thus summation 
in Vi does not appear in eq. (31). On the other hand, in eqs. (32) and (33), summation 
with respect to r j is necessary, since Ci and C 2 need not reside on the site r, where x-ray 
absorption and emission occur. All of the above formulae are valid also for general cases 
including several transition-metal atoms. 

Here we should note that the total RIXS intensity W (g, e; g', e') does not equal the 
sum of the resolved intensities, e.g., W{q, e; g', e') 7 ^ W^^^q, e; g', e')-|-hF(^’^(g, e; g', e') + 
hF('^)(g,e;g',e'), W^\q, e; q', e') ^ lT^^f\^^(g, e; g', e'), lT('^)(g, e; g', e') ^ 

and so on. This is because the total summed spectrum 
hF(g, e;g', e) contains interference terms such as fci; g, g')F‘^'^^*(rj; g, g'), while 

resolved spectra W^^\q, e; q', e') and hF^'^^(g, e; g', e') contain only (rp fep g, g')p 
and |F‘^'^^(rj; g, g')p, respectively. Nevertheless these resolved spectra are useful to take 
microscopic insights into the mechanism of RIXS, as we shall see later. 

For numerical calculation of eq. (19), we use the Lorentzian expression for the 6- 
function: 


5(12 -|- Eaj^{ki) — Ea2{ki + Q)) 


(34) 


7r[12 + E,,(fci)-E„,(fci + Q)]2 + e2- 
This function possesses poles at 12 = Ea 2 {ki + Q) — Ea^{ki)±ie, which correspond to the 
transition from band Oi to band 02 . Therefore, at a hrst glance, one might consider that 
eq. (19) describes only simple band-to-band transitions and fails to describe local d-d 
transitions. This naive view is not correct, as explained next. We should note that, for 
overall consistency, the factor e should equal the damping rate of the excited electron- 
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hole pair, already introduced above. The position of the pole = Ea^iki + Q) — 
Eaj^{ki) ± ie is modihed to a non-trivial position by the RPA correction. The modihed 
poles describe bound states between the excited electron and hole in the hnal state. As 
a result, we shall see not only charge-transfer excitations but also magnon excitation 
and local d-d excitations can be described within our HF-RPA calculation on the basis 
of electronic bands. 

For later discussions, we dehne the following quantity. 






oc 


M=-J 

.1 


{d^\e ■ r\2pjM){2pjM\e' ■ r\dp) 


(35) 


M=-J 

where ( = {i, a), (' = (£', a'), and £ and £' lie on transition-metal site r,. We should note 
that ( [C'] represents the orbital and spin of the d electron which is initially promoted 
from a 2p state [hnally decays down to the empty 2p state], and this electron could 
be different in general from the electron left above Ep [removed below Ep] in the hnal 
state at f = -|-oo. Nevertheless, we shall see e; qr', e') is useful for understanding 

the dependence of RIXS spectra on the scattering geometry. 


3. Numerical Results 

3.1 Overall spectral structure 

Throughout our present study, we set the incident photon energy at the Cu-L^ edge 
(see Appendix B for absorption spectra), and set the azimuth angles as 0 = 0' = 0. 
Firstly, we present a typical calculated result of RIXS spectrum within our theoretical 
framework (see the thick solid curve in Fig. 3). The spectrum consists of three main 
features: magnon excitation at low energies up to 400 meV, d-d orbital excitations 
around 1.7-2.1 eV, and the charge-transfer excitations between 3-9 eV, which are semi- 
quantitatively consistent with experimental spectra. 

The origin of each feature is clarihed by using the process-resolved spectra (see § 2.2, 
and eqs. (28), (29) and (30)). The low-energy d-d and magnon excitation weights orig¬ 
inate solely from the d-process. We should note that, in our electronic band approach, 
the d-d excitations are described as bound states of electron and hole on the bands. 
Therefore d-d excitations are not contributed from the Oth-order process but from the 
d-process. The Oth-order process contributes only above 2.2 eV excitation energy to the 
RIXS weight, since the insulating gap is 2.2 eV. On the other hand, high-energy (above 
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Energy Loss (eV) 


Fig. 3. (Color online) Typical calculated spectra as a function of energy loss. In-plane momentum 
transfer of photons is set to Q = q — q' = (7r/2,0). Bragg angles are 9 = 9' = 7r/4 (specular geometry). 
Initial incoming photons are in tr polarization (i.e., = 7r/2) and the final outgoing photons are in 
TT polarization (i.e., Vj' = 0). The thick solid curve is the total spectra. The thin dashed, dotted and 
solid lines are the process-resolved spectra, (q, e; q', e') (Oth-order process), VF(P^(q, e; q', e') (p- 
process), and IF^'^)(q, e; q', e') (d-process), respectively. W^^\q,e\q',e') is negligibly small. The thick 
solid and thin solid curves almost completely merge at the low-energy region below 2.2 eV. 


3 eV) part is contributed from both the Oth-order process and d-process. The p-process, 
in which the d electrons are excited by screening the inner-shell 2p core hole, is negligi¬ 
ble in the present case. As intuitively understood by the local picture, the Cu-d shell is 
completely hlled with 10 electrons in the intermediate state, and the core-hole cannot 
be screened by the Cu-d electrons. In this point, L-edge RIXS substantially differs from 
A-edge RIXS. In A-edge RIXS, screening of the inner-shell Is core hole is the essential 
process.Thus the Slater-Condon integrals between the 2p and 3d states, A°(p, d) and 
F‘^{p,d), are not signihcant parameters in the present Cu L-edge case. For the same 
reason, G^{p,d) and G^{p,d) are not signihcant. 

3.2 Magnon excitation 

Here we focus on the low-energy feature below 0.5 eV. As we have mentioned, 
magnon (spin-wave) mode appears below about 350 meV, and is dispersive as shown in 
Fig. 4(a). The theoretical dispersion of the magnon is quite consistent with experimental 
one in quantitative level as shown in Fig. 4(b). This magnon mode originates from spin- 
hip processes in the d^ 2 _y 2 orbital. Here note that we have started not from a Heisenberg 
spin model but from an electron band model to derive the magnon spectrum. To ht 
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(a.) RIXS Intensity (arb. units) 



(b) 



Q 


Fig. 4. (Color online) (a) Low-energy part of calculated RIXS spectra is depicted along the sym¬ 
metry lines in the Brillouin zone, (b) Gray-level map of the low-energy RIXS intensity. The scattering 
geometry for calculation \s 0 = 9' = 7r/4 and tt cr' (■!/) = 0, The plots are experimental 

data of magnon excitation energy, which are read from Refs. 10 and 14. Experimental data contain 
those oiO -\- O' = 90° (for small |(5|) and 130° (for large |(5|). 


well to the experimental magnon dispersion, Heisenberg-model approaches require not 
only nearest-neighbor exchange J but also long-range ones. The exchange parameters 
are often used as htting parameters. On the other hand, in our band approach, effects 
of the long-range magnetic correlations are taken into account through the inter-site 
hopping parameters, which are determined by the hrst-principles band calculation. It 
is notable that the magnon intensity should increase divergently towards <5 = 0. This 
divergence is not a result from ordinary elastic processes, since we have already excluded 
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Fig. 5. (Color online) Polarization dependence of the low-energy part of calculated RIXS spectrum. 
9 = 9' = 7r/4, and in-plane momentum transfer of photons is set to Q = q — q' = (7r/2,0). Initial 
incoming photons are in tt polarization (i.e., "0 = 0) in (a) and in a polarization (i.e., ^ = 7r/2) in (b). 


the elastic contributions, which are expressed by diagrams whose upper normally-time- 
ordered and lower reversely-time-ordered branches are disconnected (See Fig. 2(1)). 

Photon-polarization dependence of the magnon excitation is so drastic, as shown in 
Fig. 5. Intensity of the magnon excitation is maximized, when polarization directions 
of incoming and outgoing photons are perpendicular to each other (tt —)■ a' or a —?■ tt'). 
On the other hand, the weight of the magnon excitation vanishes, when polarization 
directions of incoming and outgoing photons are identical to each other (vr —)■ tt' or 
a a'). Thus, it may become possible to identify the magnon weight more clearly, 
by sorting out not only incoming photons but also outgoing photons with respect to 
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polarization. 

3.3 d-d excitations 

Here we focus on the d-d excitations around 2 eV energy loss. In Fig. 6, a calculated 
spectrum for a typical scattering geometry is compared with the experimental data 
read from Ref. 34. Two-peak structure is well reproduced by the present calculation, 
although the peak positions depend on the modihcation of the completely-hlled Cu-d 
levels (see Appendix A). In that scattering geometry, the low-energy peak around 1.7 
eV are attributed mainly to the d-d transition from d^y to dx 2 _y 2 , while the high-energy 
peak around 2.2 eV is to the d-d transition from dy^ to dx 2 _y 2 (Note that we are working 
with the electron representation not with the hole one). The transition from d^z to dx‘i_y 2 
is absent in that scattering geometry. The interval between the two peaks is somewhat 
overestimated. We consider that this overestimation is ascribed to the band structure 
calculation (and subsequent Wannier htting and Hartree-Fock calculation), since we 
retained the energy level splitting between the d^y and dy^/^z levels in modifying the 
completely hlled Cu-d levels (see Appendix A). In Fig. 6, we should note that the peaks 
are not contributed only from spin-conserved excitations but also from spin-flipped 
excitations. In particular, the spin-conserved and spin-flipped transitions from dxy to 
dx^-y^ almost evenly occur. Energy positions of d-d excitation peaks do not signihcantly 
depend on momentum transfer, reflecting their localized nature (Momentum-transfer 
dependence is not shown). 

The d-d excitation spectrum signihcantly depends on the scattering geometry. In 
the right panels of Fig. 7, the d-d excitation spectra are displayed for various Bragg 
angles with 9 -\- 6' = 7r/2. For vr —)■ a', d-d excitation intensity monotonically increases 
with increasing 9 (left panel of Fig. 7(a)). For a tt', intensity of the 1.7 eV [2.2 
eV] peak decreases [increases] drastically with increasing 9 (left panel of Fig. 7(b)). 
These behaviors seem well correlated with e] q', e')\ As seen in Fig. 6, 

the 1.7 eV peak originates mainly from the d-d transition xy ^ and the 2.2 

eV peak originates from the d-d transition yz ^ x^ — y"^. For tt a', P^^^y 2 „)(^xy „) 
and P(^^_y 2 -cr) both increase with increasing 9 (right panel of Fig. 7(a)). For a —)■ 
p, P(x^-y2a)ixya) decreases and P^^^y 2 ^)(^y^^) increases with increasing 9 (right panel 
of Fig. 7(b)). In the left panels of Fig. 8, the d-d excitation spectra are displayed 
for various polarizations of the incoming and outgoing photons. The dependence on 
polarization can roughly be understood in the same manner as above. For example. 
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Fig. 6. (Color online) A typical d-d excitation spectrum around 2 eV. Thick and thin solid curves 
denote the calculated results of total and spin-orbital-resolved RIXS spectra, respectively. For the 
calculated results, polarization of photons is set as tp = tp' = ttI2 (tr — >■ tr'). For the weight from xy —>■ 
x^ — y^, four curves {xy t— >■ x^ — y^ ti x^ — y^ I, xy J.—>■ x^ — y^ t; xy x^ —y'^ ),) almost evenly 

contribute and their curves are overlaid (Note that we are working with the electron representation 
not with the hole one). Circles are experimental data read from Ref. 34 (incident incoming photons are 
in a polarization, but polarization of final outgoing photons is not discriminated in the experiment). 


when the incident photon is in tt polarization ('0 = 0), intensity of the 1.7 eV [2.2 eV] 
peak increases [decreases] drastically with increasing ip' (left panel of Fig. 8(a)). This 
behavior is understood well by noticing that ^ increases and ^ 

decreases with increasing ip' (right panel of Fig. 8(a)). 

4. Discussions 

In this section, we present some remarks on L-edge RIXS, our formulation and 
calculated results. 

As we have seen, in L-edge RIXS, the d-d excitations play a major role in the RIXS 
weight. This is in strong contrast to the case of iC-edge RIXS. In the case of iC-edge 
RIXS, the process of screening the created Is hole by the d electrons is dominant, where 
the d electrons are scattered isotropically due to the strongly localized Is core hole and 
usually prohibited to change their orbital state. In those processes, the d electrons 
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(a) \|f=0, \|;’=Ji/2 (Ji-^o’) 




(b) \|f=3i/2, \|f’=0 




Fig. 7. (Color online) Dependence of the d-d excitation spectrum (left panels) and 

(I'ight panels) on the Bragg angle 6 (or dependence on d = {0 — O') 12) 
with maintaining 9 -\-O' = 7r/2. 0 = tt/S, 7 r/4 , Stt/S correspond to 5 = —tt/S, 0,+7 r/8, respectively. 
In the right panels, e] q', e') is shown for P = xy,yz,xz,iz'^ — . Empty and filled 

symbols are for the spin-conserved (cr = u') and spin-flipped (ct = —a') components, respectively. In 
(a), photon polarizations are set as '0 = 0 and if)' = 'x 12 (tt —s- a'). In (b), photon polarizations are set 
as 0 = x/2 and -0' = 0 (cr —?► tt'). Inset in (a) schematically represents the scattering geometry for the 
three incident Bragg angles, 0 = tt/S, 7 r/4, Stt/S. 


conserve their orbital angular momentum. This means that the orbital-diagonal d-d 
excitations are dominant in iC-edge RIXS. As a result, completely hlled d orbitals in 
the initial state are completely hlled still in the hnal state of RIXS, and completely 
empty d orbitals are empty still in the hnal state. Therefore we can construct a good 
ehective theory by neglecting completely hlled or empty d orbitals, when we describe K- 


/ »-» -I 
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(a) \|;=0 



(b) V|;=ji/2 




Fig. 8. (Color online) Dependence of the d-d excitation spectrum (left panels) and 
-F’(a;2_'j,2 (right panels) on polarizations of incoming and outgoing photons. In the 

right panels, cr)(£'cr') (?’ is shown for P = xy, yz, xz, 3z^ — Empty and filled symbols 

are for the spin-conserved (a = a') and spin-flipped (a = —<j') components, respectively. In (a) [in 
(b)], the incident incoming photons are in tt [cr] polarization, i.e., tjj = 0 [■(/> = 7r/2]. The Bragg angles 
are 9 = 9' = tt/A in the both cases. 


edge RIXS. In fact, orbital excitations between different d orbitals seem inactive in most 
cases of /C-edge RIXS, and we can take only active d orbitals, e.g., the dx 2 _y 2 orbital 
in the case of cuprates with the tetragonal or octahedral coordination. On the other 
hand, in L-edge RIXS, orbital states of d electrons can change in the hnal state. This is 
because the d electrons can exchange their orbital angular momentum into spin angular 
momentum and vice versa, mediated by the 2p states in which spin and orbital angular 
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momenta are coupled. Therefore, to analyze L-edge RIXS spectra over a wide-energy 
range, we should include d orbitals which are initially completely hlled or completely 
empty, in addition to partially hlled d orbitals. In this sense, L-edge RIXS is essentially 
a multi-orbital phenomenon, and can involve much more complex excitation processes 
than iL-edge RIXS. 

The p-process, where the 2p core hole is screened by the Cu-d electrons, is almost 
ineffective in the present case of Cu L-edge RIXS. This is because the narrow Cu-d 
band is completely hlled in the intermediate state of RIXS, and is in great contrast to 
the case of L'-edge RIXS. However, we should note that the p-process can be ehective 
in other transition-metal compounds whose transition-metal d bands are broad (i.e., 
localization is weak) and are not completely hlled in the intermediate state. 

As we have seen in § 3.2, the magnon intensity becomes divergent towards Q = 0. 
Very recently, Igarashi and Nagao found that such divergence essentially appears around 
<5 = 0 as a result from broken symmetry.^®) According to them, this divergence orig¬ 
inates from anisotropic terms of the scattering amplitude which are a natural result 
from the long-range ordering, and is not reproduced within the so-called fast collision 
approximation (FCA). Although it is beyond the scope of the present work and not 
straightforward to clarify the analytic relation between our present calculation and 
theirs, such anisotropic terms are included within the HF approximation in our present 
calculation. Our calculation suggests that this divergent behavior should survive for 
realistic values of core-hole life time. Such divergent behavior around Q = 0 is a sub¬ 
stantial difference between RIXS and neutron scattering, and deserves to be searched 
for carefully, although the measurement should be performed at low temperatures well 
below the Neel temperature and the inelastic weight around Q = 0 may be difficult to 
distinguish from the elastic one experimentally. 

Note that the multi-magnon processes are beyond the present theoretical approach, 
n-magnon generation processes are represented by diagrams in which the upper and 
lower branches are connected by 2n lines, which are omitted in Fig. 2(1). Inclusion of 
the multi-magnon excitation processes is not feasible for the present complex electronic 
structure and remains an interesting future work. 

In § 3.3, we have interpreted the dependence of the d-d excitation spectrum on 
the scattering geometry by using the square of the product of the dipole-transition 
matrix, e; g', e') (see eq. (35)). We should note that such interpretation is not 

always valid, since RIXS intensity is not determined only by the dipole-transition matrix 
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elements. Analysis based on e; q', e') will become less effective for systems where 

orbitals are strongly hybridized with each other due to low crystal symmetry. In general, 
the d electron which is initially promoted from a 2p state can be quite different from 
the d electron remaining above Ep in the final state, as a result from the intermediate 
scattering processes and orbital hybridization. 

5. Conclusion 

We have discussed the Cu L^-edge RIXS for a typical parent compound of high-Tc 
cuprate superconductors La 2 Cu 04 on the basis of the first-principles electronic bands. 
We consider that the advantages of our approach are: it can reproduce RIXS features 
in the wide-energy range, including the low-energy magnon spectrum, intermediate- 
energy d-d excitations and high-energy charge-transfer excitations, consistently with 
experiments. It is applicable to relatively complex three-dimensional multi-orbital sys¬ 
tems. Finite size effects are absent. Therefore it is useful for analyzing experimental data 
for a wide variety of transition-metal compounds. Disadvantages are that the ground 
state needs to be described accurately within the Hartree-Fock approximation, and 
multi-magnon weights are not included at all. 

In conclusion, we would like to stress the importance of photon polarization depen¬ 
dence: the intensity of magnon excitation and the spectral structure of d-d excitations 
depend significantly not only on the polarization direction of incident incoming photons 
but also that of outgoing photons. Polarization of scattered photons is not discriminated 
so far in most of experiments.^®^ A full knowledge of the both polarizations of incoming 
and outgoing photons could become crucial to correctly assign RIXS spectral weights 
to magnetic and orbital excitation processes. 
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Appendix A: Hartree-Fock Approximation 

By fitting to the first-principles electronic structure of the nonmagnetic state, we 
obtain transfer integrals and one-particle energy levels. We take about 2000 transfer 
integrals to fit 17 bands near the Fermi energy using the 17 Wannier orbitals (all of the 
Cu-d and 0-p states are included). Some typical values of calculated transfer integrals 
are 1.346 eV between nearest-neighbor Cu-dx 2 -y 2 and in-plane O-p^/y, and 0.623 eV 
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between nearest-neighbor in-plane O-pxjy These values quantitatively agree with 1.3 
eV and 0.65 eV in Ref. 37. One-particle energy levels are obtained for the Cu-d orbitals 
as: Exy = 6.81 eV, Syz/xz = 7.12 eV, ex 2 -y 2 = 7.55 eV, = 7.56 eV (the Fermi 

level is about 9.53 eV). Only the dx 2 -y 2 state is partially hlled, and the other d orbitals 
are completely hlled with electrons. Here we consider that these values do not rehect 
a realistic level scheme of the local Cu-d orbitals: for example, Sx^-y^ — Syzjxz = 0.43 
eV will be too small. To reproduce the d-d transition energies which were recently 
determined with great accuracy by L-edge RIXS,^^^ we modify the one-particle energy 
levels of the completely-hlled Cu-d orbitals in the following way: we subtract 1.40 eV 
from Si for the t 2 g {de) orbitals and 1.65 eV for the d^^z^-r^ orbital, and as a result we 
obtain Sxy = 5.41 eV, Syz/xz = 5.72 eV, Sx^-y^ = 7.55 eV, e^z^-r^ = 5.91 eV (the Fermi 
level is about 9.49 eV). Hereafter, we redehne ei by the subtracted one-particle energy. 
In Fig. A-l(a), we show the original electronic structure, the htted bands using the 17 
Wannier orbitals, and the bands obtained for the modihed Cu-d energy levels. 

To describe the antiferromagnetic ground state, we apply the HF approximation to 
the tight-binding Hamiltonian 77^./.- For fhe Coulomb integrals introduce 

the following notation: 


(A-l) 

(A-2) 



Ji^i! and are the so-called direct and exchange integrals, respectively. We assume 
spin polarization is induced only in the d orbitals at each Cu site, and take mean helds 
only for the Cu-d electrons. The mean-held Hamiltonian for H^j. is 





i I i i 


L.lll. j L.lll. 'idili jy / \ 

-XE 2 {nii){niii) + EE f {nii){nui) -f {rriii) ■ (mj^/) j , (A-3) 



where 



(A-4) 


a 
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mu = (A-5) 

crcr' 

using the Pauli matrix vector <t. Within the HF theory, we should consider that the 
one-particle energy si is already including the following energy shift from the bare one, 

^hf ^ ^ (A-e) 

due to the electron-electron Coulomb interaction at transition-metal site rj. Therefore, 
before determining the magnetic ground state, we need to evaluate the bare one-particle 
energy by = ee — Af^, where Af^ is evaluated from the expectation values of parti¬ 
cle numbers {nii)^s in the paramagnetic state using eq. (A-6). Note that Sis here are the 
modihed values mentioned above. For the Coulomb integrals given in § 2, the obtained 
values of are as follows: e^xy = —39.69 eV, = —40.02 eV, £^^ 2 _y 2 = —38.86 

eV, £^^^ 2 _r 2 = —40.06 eV. Maintaining these values of we determine the mean-helds 
{riii) and (mu) self-consistently. We assume the antiferromagnetic ordering with the 
spin moments pointing along [110], as observed in neutron scattering.As a result, 
we obtain 68 diagonalized energy bands {Ea{k), 1 < a < 68) for the antiferromagnetic 
ground state (Note spin degeneracy and folding of the Brillouin zone). The dx^-y^ or¬ 
bitals accommodate 1.27 electrons and take the spin moment of 0.69 /i^ per Cu site, 
while the other Cu-d orbitals are completely hlled with electrons. Figure A-l(b) shows 
the obtained density of states for the antiferromagnetic ground state. 


Appendix B: Resonant X-ray Absorption Spectrum 

Resonant x-ray absorption spectrum (XAS) is obtained by calculating the number 
of photons absorbed by the electronic system per unit time. Keldysh diagrammatic 
representation of XAS is presented in Fig. B-l(I), and the analytic expression is 




t.m.u. @ri 


dt e 


—iujt 


EXE E 


i CC' J M=-J 

WcjMiru q, q, 6 )^++ q, -t), (B-1) 


where q = (cn, q) and Keldysh Green’s functions 

for the 2pjM electrons at transition-metal site i and that for the transition-metal d 
electrons, respectively: 
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(a) 




Energy (eV) 


Fig. A'l. (Color online) (a) Electronic band structure of the non-magnetic state for La 2 Cu 04 . Thick 
solid, thin solid, and dashed curves are the results of the band calculation, Wannier-function fitting, 
and the tight-binding model with modified Cu-d levels (see the text for details), respectively. The 
Wannier-htting curves are almost completely overlaid on the 17 curves of the band calculation result 
around the Fermi energy, (b) Density of states (DOS) for the antiferromagnetic ground state. Thick 
solid, thin dashed, and thin solid curves are the results of the total, t 2 g and Cg partial DOS, respectively. 
The Fermi energy is set to zero in the both panels. 


a 

Thus eq. (B-1) is reduced to 

t.m.u. @ri J 

A{q,e) = 2ti 

i CC' J M=-J 

+ e2pj{ri) - Ea{k)). (B-4) 

k a 

































J. Phys. Soc. Jpn. 


(I) 


( 11 ) 



Fig. B-1. (I) Diagrammatic expression of x-ray absorption intensity. Wavy line, thick solid line, and 

thin solid line represent the propagators for the photon, Cu-2p and Cu-3d electrons, respectively. The 
shaded triangle is the electron-photon interaction renormalized by the Coulomb interaction between 
the 2p and Sd electrons. (II) The vertex correction by the multiple scattering between the 2p and Sd 
electrons (Ladder approximation). 




q, e) is the effective 2p-d dipole-transition matrix renormalized within the 
ladder approximation as in Fig. B-l(II): 

@ri J' 

w^,jM{ri]q,e) = q, e)-b z ^ ^ ^ V2p-d{ri] J'M'X'X", JM) 

C'C" J' M'=-J' 

1 r+cxD 

e) — ^ + q,z + u) 

k ^ 

@ri J' 

= wc,JM{r^■,q,e)-J2J2 E ^2p-d(n; J'M', C; C", ^M) 

C'C" J' M'=-J' 


X 




1 - na{k) 


k a 


UJ+ £2pjX'f'i) - Ea{k) +ie 


(B-5) 


where G 2 pjj^(i){k) and (j^/^^(/c) are the ordinary causal Green’s functions for the 2pjM 
electrons at transition-metal site i and that for the transition-metal d electrons, respec¬ 
tively. 

Calculated results of XAS for several values of the Coulomb interaction between 
the Cu-2p and Cu-Sd orbitals are displayed in Fig. B-2. The main peak position shifts 
to the low-energy region, increasing the Coulomb interaction. The calculated spectral 
shape suggests a possibility that a small satellite peak can appear. For the calculations 
of RIXS spectra, we simply use the unrenormalized matrix elements tPc,jM(^i; q, e). 
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Fig. B'2. (Color online) Calculated x-ray absorption spectra (XAS) near the Cu ia-edge. We show 
four cases of the Coulomb interaction between the Cu-2p and Cu-3d orbitals: d) = d) = 

0,1, 2,4 eV. The Bragg angle is 6 = 7r/4, and the absorbed photon is in tt polarization {ip = 0). 


The photon energy for the main resonance absorption is changed by using renormalized 
Q, e), but it can be re-adjusted to the observed value by shifting the inner-shell 

2p level. 
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